o 

(N 
(N 



On the lateral migration of a slightly 
deformed bubble rising near a vertical plane 

wall 

By KAZUYASU SUGIYAMA^ and FUMIO TAKEMURA^ ^ 

^Department of Mechanical Engineering, School of Engineering, The University of Tokyo, 

7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan 

^National Institute of Advanced Industrial Science and Technology, 1-2-1 Namiki, Tsukuba, 

Ibaraki 305-8564, Japan 



O 



oo 



Deformation-induced lateral migration of a bubble slowly rising near a vertical plane 
wall in a stagnant liquid is numerically and theoretically investigated. In particular, our 
focus is set on a situation with a short clearance c between the bubble interface and 
the wall. Motivated by the fact that numerically and experimentally measured migration 
velocities are considerably higher than the velocity estimated by the available analytical 
solution using the Faxen mirror image technique for a/{a + c) ^ 1 (here a is the bubble 
radius), when the clearance parameter e{= c/a) is comparable to or smaller than unity, 
the numerical analysis based on the boundary-fitted finite-difference approach solving 
the Stokes equation is performed to complement the experiment. The migration veloc- 
ity is found to be more affected by the high-order deformation modes with decreasing 
e. The numerical simulations are compared with a theoretical migration velocity ob- 
tained from a lubrication study of a nearly spherical drop, which describes the role of 
the squeezing flow within the bubble-wall gap. The numerical and lubrication analyses 
consistently demonstrate that when e < 1, the lubrication effect makes the migration 
velocity asymptotically ^Vgi/{25e^) (here, Vbi, /i, and 7 denote the rising velocity, the 
dynamic viscosity of liquid, and the surface tension, respectively). 
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1. Introduction 



Recent technical progress in generating microbubbles (e.g. iGarstecki et al\ ( 20061 ): 
Makuta et all ( 20061 )). including potentials as actuator and sen sor, has enhan c ed th e 



range of applications, e. g. additives to reduce a turb ulence friction (ISerizawa et al. (12005 )') 



drug d elivery capsules ( Shortencarier et al\ ( 20041 )). and contrast agents ( Correas et al 
(|2001l )). In many situations, a bubble encounters a boundary wall during its transport 
process, and a hydrodynamic interaction occurs as characterized by the inter-scalc be- 
tween the bubble and the wall. In practice, it is of primary importance that the bubble un- 
dergoes a r epuls i ve or a t tractive force in the wall - normal direction, whi ch causes a lateral 
migration (JLeall (|l980h : iMagnaudet et al\ (|2003[ ): iHibiki fc Ishi3 (|2007[ )) and determines 
the bubble distribution, when translating parallel to the wall. As the simplest model sys- 
tem, one may raise a phenomenon of a spherical bubble rising near a vertical infinite plane 
wall in a creeping (Stokes) flow. However, there is no mechanism to generate the lateral 
migration force, as kinematic revers ibility is en sured by symmetry of the boundary and 
by linearity in the Stokes equation (ILeall (|l992| ). chapter 4). In f act, the migration forc e 



by linearity in the Stokes equation (ILeall (|l992| ). chapter 4). In f act, the migration fore 
stems frorn nonlinearities in the advectiye morn e ntum transport ( Cox fc Brenneij ( 1968 
iHo fc Leall (flQTJ : I Vasseur fc Coxl (|l976L Il97lt ): ICox fc Hsul (|l977l ): iMcLaughlh] (|l993 
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Figure 1. Schematic of a buoyant bubble moving near a vertical plane wall in a quiescent 
liquid. The solid outline shows a deformed interface with a deflection / from a spherical interface 
indicated by the dotted outl i ne, (a ) left panel: coordinates system around the bubble. In the 
analysis of lMagnaudet et al.\ (|2003l ). a/d is assumed to be sufficiently smaller than unity (here, 
a is the bubble radius and d is the distance between the bubble centroid and the wall), (b) right 
panel: inter-scale coordinates between the bubble surface and the wall and scaling relations 
suitable to a lubrication theory. A clearance parameter is defined as e = c/a. 



Cherukat fc McLaughliJ (Il994: ISecker et al\ (Il996l): iMagnaudet et al\ (120031)) and/or 
the interfacial deformabilitv (IChafFev et al. { 1965 k Chan fc Leal ( 19791): Shapira fc Habei 



lie interracial aerormapiiitv liunattev et ai.i iii^boii: unan az: i^eaJllliiJ/aKlonapira &; Mapeil 
1988|) : IUiitewaal et ad(|l993t ): IUiitewaal fcNiihoff(|l995l ): IMagnaudet e^ ad(|2003[) : IWang 



(120061)), to break the symmetry. 



For a tank-treading vesic l e tran slating parallel to the wall, the migration force was the- 
oretically obtained by lOllal ( 1997 ). in which the vesicle shape is prescribed as a strongly 
non-sp herical ellipsoid, and the theory was experimentally validated by ICallens et al. 
( 20081 ). For a bubble or drop, the shape cannot be prescribed since it obeys the Laplace 



law and depends on the surrounding fluid flow. The theoretical success in solving the non- 
trivial problem of the de formation-induced migration of the bubble or drop was made by 
Magnaudet et al\ (J2003l) using the Faxen mirror image tec hnique and the Lorcntz recip- 
rocal theorem. However. I Wang fc Dimitrakopouloa ( 20061 ) performed a numerical study 
on the motion of a drop with the same viscosity as the surrounding fluid in a linear 
shear flow by means of a boundary element method, and pointed out that the theory 
considerably underestimates the migration velocity or erroneously predicts the lateral 
motion, despite c onsistent predictions o f the rising velocity and the interfacial defor- 
mation. Recently. iTakemura et al\ ( 20091 ) experimentally measured the lateral migration 
velocity of slightly deformed bubbles in a wall-bounded shear flow, and found a clear 
discrepancy between the experimental and theoretical values of the deformation-induced 
transverse force. Then, they computed the quasi-steady evolution of deformable bubbles 
moving in a wall-bounded linear shear fl ow at zero Reyno l ds nu mber using a spectral 
boundary clement method developed by iDimitrakopoulosI ( 2007 ). and found that the 
measured deformation-induced lift force agrees quantitatively well with the computa- 
tional prediction. Motivated by their conclusions, we revisited exp e rimen tal data of the 
bubble migration in a quiescent liquid obtained by ITakemura et al\ (2002), and analysed 
the data on the conditions that the clearance between the bubble interface and the wall 
is comparable to or shorter than the bubble radius, which was not considered there. The 
results revealed that the discrepancy between the migration velocities of the experiment 
and the theory increased as the bubble was closed to the wall, as detailed below. 

In this paper, we focus on the bubble motion in a quiescent liquid to simplify the 
subject. Let us consider the migration velocity Vb3{= Vb • 63) of a bubble rising near a 



Dimitrakopoulos 



On the lateral migration of a slightly deformed bubble 



0.6 ■ 



u 

^ 0.2 



0.0 




o; 

Bo=0.31 
(0.066<Ca<0.083) 

\/ o/ 



o/ 



0.0 0.2 0.4 0.6 0.8 1.0 

Bo=0.19 
(0.046<Ca<0.055) .->' 



■ o / 



-e^^ 




0.0 



0.2 



0.4 



0.6 



0.8 



K- 



1.0 



Figure 2. Migration veloci t y Vb3 versus k{ 

solution (|Magnaudet et al\ (|2003l ')'l VBs/iCaVBi) = 3k^ 



(l + £) 



^). The solid curve indicates the analytical 
" (1 + 3k/2) /{40(1 + 3k/4)} with the 
assumption of the sufficiently long distance between th e bubble centroid and th e wall, i.e. e 2> 1. 
The dashed curves indicate the experimental results (|Takemura et al\ (|2002f )). the circles the 
results obtained by the numerical simul ation, and the dashed-dotted curve the prediction (|2.22p 
by means of the lubrication approach (|Hodges et al.\ (|2004l )'). The inset shows the ratio of the 
simulation result to the analytical solution. 



vertical plane wall at a distance d between the bubble centroid and the wall as schemat- 
ically illustrated in figure [Ija). The bubble has an equivalent radius a to that of a 
sphere with the same volume. Introducing an interfacial deflection f{9, (j)) from a sphere, 
we write the distance from the bubble centroid to the interface as a + f. The experi- 
mental results used here were measured under the condition that the Reynolds number 
Re = 2paVBi//i (here Vbi (= Vb • ei), p, and ^ respectively denote the rising veloci ty, the 
density, and the dynamic viscosity of liquid) is unity or less ( Takemura et ali ( 20021 )). The 
pure lateral migration velocities induced by the deformation Vb3 were calculated from 
the measured values substituting the velocities induced by the inertia effects. Follow ing a 
Stokes flow theory for the deformation- induced migration ( Magnaudet et all ( 2003[ )). we 
can characterize the system using two parameters, i.e., a clearance parameter e{= c/a) 
and a capillary number Ca = ^Vbi/i (or a Bond number Bo = pa^g/^ as used in 



Magnaudet et al. 



( 20031 1). Here 7 and g respectively denote the surface tension and the 
acceleration of gravity. Further, as long as Ca <C 1, we may use Ca as a perturbation 
parameter, and reduce the Ca-dependent system to another, in which Vbi, Vsa/Ca and 
//Ca are dependent only upon e, under the infinitesimal deformation assumption. Figure 
[2] shows the migration velocity Vb3 away from the wall normalized by CaVei as a func- 
tion of k{= (1 -I- e)^^). (It should be noticed that although k as well as e are measures 
of the distance between the wall and the bubble, hereafter k is also used to make some 
equations for the wide gap case {k <C 1) simple.) The measured velocity is found to be 
much higher than the analytical solution especially for the large k. A possible inference 
drawn from this result is that there exists an additional ingredient to generate repulsive 



force f or narrow bubble-wall gap, which is not covered by the theory of lMagnaudet et al. 
( 200,lh . 



As the most crucial restriction involved in the mirror image technique, we can raise an 
assumption that the bubble- wall distance is much longer than the bubble radius. However, 
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with regard to an inertia effect on the lateral velocity of a rigid sphere, iTakemural ( 20041 1 



experimentally demonstrated that the mirror image approach has robust applicability 
in prediction beyond the wide-gap precondition. At this moment, we cannot conclude 
whether the discrepancy in the deformation-induced migration velocity comes from the 
erroneous prediction due to the miscalculation or the contradictory conditions in the 
boundary element computation and the experiment with the theoretical assumptions. 
To complement the experiment and to gain further insight into the deformation-induced 
migration, we investigate the m igration behavior with at tention to shortness of a bubble- 



wall clearance c(= d — a). As in lMagnaudet et al\ (|2003f ). using two Stokes flow solutions 



for a spherical bubble translating parallel and perpendicular to the wall, we apply the 
Lorentz reciprocal theorem to evaluating the migration velocity. We carry out numerical 
simulations using a boundary-fitted grid, which can accurately implement the boundary 
conditions and release the constraint of the sufficiently wide bubble- wall gap in the mirror 
image technique. In addition, comparisons with a theoretical migration velocity for e <S^ 1 
obtained from a lubr ication study of a nearly spherical drop moving near a tilted plane 
( Hodges et all ( 20041 ) ). in which the secondary flow due to the change in the boundary 



geometry caused by the bubble deformation is responsible for the wall-normal force, are 
made to shed more light on the short clearance effect. 



2. Numerical simulation 

2.1. General formulation 
To clarify the physical mechanism of the repulsive force , we nu merically address the bub- 



ble migration. In a similar manner to lMagnaudet et al\ (J2003I ). instead of directly solving 



the flow field with the deformed bubble, we employ the Lorentz reciprocal theorem to de- 
termine the lateral migration force and velocity through coupling two fiow fields around 
a spherical bubble translating parallel and perpendicular to the wall. In the subsequent 
developments, the basic equations and the involved variables are nondimensionalized us- 
ing a, Vbi and fi. We suppose that the bubble quasi-steadily rises near an infinite fiat 
plate in a stagnant incompressible liquid, and both the Reynolds and capillary numbers 
are sufficiently smaller than unity. Hence, the system is described by the steady Stokes 
equation for solenoidal velocity vectors, i.e. 

V • C/ = V ■ M = 0, —CO < Xi < OO, — OO < X2 < oo, — 1 — £ < ^3 < oo, (2.1) 

V • S = V • (T = 0, — oo < xi < oo, — oo < X2 < oo, — 1 — e < xs < oo, (2-2) 

where (C/, S) and (it, a) are the velocity and stress fields for the bubble translating 
respectively parallel and perpendicular at a speed of unity to the wall. The Ca depen- 
dence of the interfacial deflection is given by f {9 , (/}; Ca.) = Caf^'^^\9,4>)- The bubble 
deformation obeys the Laplace law for the infinitesimal defiection |/| ^ 1 with Ca ^ 1, 



(V^ + 2) /(^^) = -n ■ S • n + 3a-i {xm ■ E • n)ss at ^ xl + xl + x^ =: 1, (2.3) 

where n represents the normal unit vector pointing outwards the liquid, Vs(= V — n{n ■ 
V)) is the nabla operator along the tangential directions on the bubble surface, {■■■)sb is 
the area average taken over the bubble surface, and xi is the coordinates in the upward 
direction from the origin at the bubble ccntroid. Kinematic and free-slip conditions are 
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imposed on the bubble surface, i.e. 

nU =0, 

rz : r 



{n-Yl)-><n =0, 
n- u =0, 

(n • cr) X n =0, 



at 



(2.4) 



at xa = —1 — e. 



where we take the reference frames viewed from the bubble. On the plane wall, we impose 
the no-slip condition 

U =-ei, 

U = 63, 

Sufficiently far from the bubble, the velocity vectors approach the uniform velocities 
U ^-ei, 

U -^ 63, 



as 



(2.5) 



(2.6) 



Thanks to the reciprocal theorem ( Leall ( 19801 ): (35) of iMagnaudet et al\ ( 20031 )). the 



(Ca) 



deformation-induced lateral force Fm = GaF^j ' to cancel the migration velocity and to 
maintain the wall-parallel motion is expressed as 



F, 



(Ca 



M 



d'x C{f 

Sb 



(Ca)N 



where Sb denotes the bubble surface, and the operator C is given by 

f dU ^. „ \ fdT, 



£ = n- cr- n ^r— n — U-\7, 
\ on 



dn 



■n-SV, 



(2.7) 



(2.8) 



— {n ■ "E ■ n — 3xi{xin ■ S • 11)3^}'^' • Vs 
The migration velocity Vbs ~ CaV^g^ is expressed as 



V, 



(Ca) 



B3 



F, 



(Ca) 



M 



Vbi 



Fdc 



where 



F, 



DC 



d-x e^ ■ a ■ n 



(2.9) 



(2.10) 



denotes the drag force acting on the bubble translating perpendicular to the plane wall. 

2.2. Simulation method 

The basic equations are numerically solved by the second-order finite-difFercncc method 
discrctizcd on the bipolar coordinates (^, 77) grid, which is boundary-fitted on both the 
bubble surface and the plane wall (see figure |9l[a) in Appendix |A|. We take care of 
the mass and momentum conservations in a discrctizcd form. For technical detail on the 
discretization, see Appendix [X] The number of grid points is N^ x Nn — 200 x 200, and the 
grid is non-uniform and refined near the wall and the bubble surface. The compu t ationa l 
procedure is based on a Simplified-Marker-And-Cell method ( Amsden fc Harlowl (1970|)) 



with a first-order Eulerian implicit time marching scheme. Such an unsteady scheme 
enables us to check whether the computation converges to the fully developed state 
through temporal changes in the budgets of the momentum and kinetic-energy transports. 
To avoid a problem associate d with singularities i n the discretization near the axis, we 
follow a method proposed by iFukagata fc Kasagil (|2002l) . 

For each run, we confirm that drag forces, numerically evaluated on the bubble sur- 
face, for both the perpendicular and parallel motions are in good agreement with the 
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respective kinetic-energy dissipation rates, numerically integrated over the entire com- 
putational domain, normalized by the translational velocities with an error of less than 
0.040%. Such an agreement between the surface and bulk quantities indicates that the 
computation is well converged and reaches to the steady state in view of the monicntum 
and kinetic-energy budgets. Further, the drag force for th e perpendicu lar motion shows 
quantitative agreement with the infinite series solution of iBarti (|l968r ) with an error of 
les s than 0.043%. Th e drag force for the parallel motion approaches the wide-gap solution 
of iMagnaudet et al\ (|2003t ) with increasing e. To make sure of numerical stability and 
accuracy, we set the clearance parameter in a range of 10^"^ < e < 9. We performed the 
convergence tests with varying the size i?max of the computational domain, the number 
iVj of nodes in the gap between the bubble and the wall, and the number Njj of nodes 
describing the bubble surface. We confirmed that the relative errors in the migration 
velocity Vb3 and the drag forces for the perpendicular and parallel bubble motions to 
those obtained on the base meshes 7V^ x iV^ = 200 x 200 for various e decrease with in- 
creasing the size i?max and the resolutions N^ and iV^. From the convergence behaviors, 
we deduce to obtain the migration velocities Vb3 with errors of much less than 1% using 
the present base meshes, which are accurate enough for the subsequent discussion. 



2.3. Migration velocity 

As shown in the inset of figure [2J the ratio of the simulation to the analytical migration 
velocity becomes close to unity as k(= (1 + e)^^) app roaches zero, indicating the simula- 
tion result is consistent with the analytical solution ( Magnaudet et al\ ( 20031 )) for small 
k{<^ 0.3). Such a consistency betwee n the different app r oache s may deny the erroneous 
prediction of the m igration velocity in Magnaudet et ali (J2003t ). of which the possibility 
was pointed out bv lWang fc Dimitrakopoulod ( 2006^ . By contrast, the simulated migra- 
tion velocity for the bubble closer to the wall with the clearance shorter than the bubble 
radius (i.e., for k >~ 0.5 presumably beyond the theoretical precondition k ^ 1) is 
considerably higher than the analytical solution. Although the simulation result reveals 
lower velocity than the experimental one, the tendency of the higher velocity than the 
theoretical one for the narrow gap is qualitatively similar to that in the experiment. Thus, 
the present simulation also indicates the presence of the additional narrow-gap repulsive 
force. 

For an undeformed spherical bubble at small but non-zero Reynolds numbers < 
Re <S 1, using the solutions to the equation set (|2.1[) - (|2.6|) . we can also evaluate the 
inertia effect on the migration velocity Vb3 



Re Vb3° from 



p(Ro) _ 



^Jd^x{es-u)-{{U-\7)U}, 



(2.11) 



V, 



(Ro) 



B3 



F, 



(Ro 



M 



(2.12) 
Vbi Fbc 

where V stands for the entire volume of liquid around the bubble. The force expression 
(|2.1ip is theoretically justifie d for the case th a t the wall is placed within a Stokes expan- 
sion region, i.e., 0(Re) < k (jVasseur fc Coxl (jl976l )). Figure [3] shows the inertia- driven 
migration velocity as a function of n. The sim ulated profile is globally consistent with the 
analytical solution (IMagnaudet et al\ ( 2003 )) even for the narrow gap k -^ 1, as opposed 
to the profile of the deformation-induced migration velocity in figure [2] It should be 
noticed that the expression (|2.7I) of the deformation-induced lateral force is written in a 
surface integral form, while the inertia-driven force p.ll|) in a volume integral form. The 
overall agreement in the inertia-driven migration velocity indicates that capturing the 
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K 

Figure 3. Inertia-driven migration velocity Vbs versus k{— (1 + £)~^) for spherical bubb le 
at < Re ^ 1. The solid curve indicates the analytical solution (jMagnaudet et al\ (|2003l )) 



FB3/(ReV'i3l) = ^(l+|K-0.516K2)(l- 



64' 



under the assumption of (0(Re) < k <C 1) 



as in figure [2] The circles indicate the results obtained by the numerical simulation 



bulk velocity distributions is important for predicting the migration velocity and can be 
robustly attained by the mirr or image technique ove r the wide range of k . The wide range 
agreement with the theories (jVasseur fc Coxl (|l976h : iMagnaudet et al\ (|2003l )) was also 
ex perimentally demo nstrated for the sedimenting rigid particle in a range of 0.1 < Re < 1 
bv iTakemural ( 20041 ) as long as the wall is placed in the Stokes expansion region. Con- 
trastingly, as shown in figure [21 the larger discrepancy between the deformation- induced 
migration velocities of the simulation and the theory with increasing k indicates that 
the migration velocity is sensitive to the local effect leading to the additional narrow-gap 
repulsive force, which may not be covered by the mirror image technique. 



2.4. Interfacial deformation 

To demonstrate the narrow- gap effect, we investigate the bubble deformation. We here 
examine the scaled interfacial deflection /(^'*) (0) = /(0, cj)) /{Ca. cos (j)). In the experiment, 
we estimated f{9, </>) taking a circumference of the bubble on the plane X2 — 0. Figure H] 
shows the angular profile of the deflection — /(^^) for various k(= (1 -|- e)~^). As shown 
in figure Ufa) for t he rel atively wide gap e = 0.67, the analytical solution for k <C 1 
(JMagnaudet et al\ (J2003[ )) is consistent with the measured and simulated deflections. 
Note that the agreement between the theoretical and simulated deflections is confirmed 
to be better in the wider separation. For the narrower gap (e = 0.10, 0.25), by contrast, 
the analytical solution of the deflection magnitude is smaller than the measured one 
especially in the wall neighborhood {6 ^ 0), which may be related to the considerable 
underestimation of the migration velocity as shown in figureO Contrastingly, the present 
simulation quantitatively captures the local near-wall profile of the measured deflection 
as well as the global magnitude. 

To quantify the local effect of such a large discrepancy in f^'^^^ on the migration 
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Figure 4. Angular profile of the deflection. The solid curves indicate the analytical solution 
IjMagnaudet et al\ (|2003n ) -/(°^) = Sr^ {i_|_|k (1+|k+|| k^)} sin 6» cos 6* with the assumption 
of the sufficiently long distance between the bub ble centroid and the wa ll, i.e., k(= (l+e)~^) <C 1. 
The circles indicate the experimental results (JTakemura et al\ (|2002l ')). the dashed curves the 
results obtained by the numerical simulation, and the dashed-dotted curves the prediction (|B 12[l 
^_j(Ca) ^ ^log(l -I- 1^)) by the lubrication approach (jHodges et al\ (|2004 l'). Upper, middle, 
and lower panels show the results at e = 0.67, e = 0.25, and e = 0.10, respectively, and the 
corresponding capillary numbers in the experiment are Ca = 0.080, Ca = 0.068, and Ca = 0.056, 
respectively. 



velocity, we here describe the deflection in an expansion form 



/(Ca)(0)^^^(Ca)p^l(^^^^ 



(2.13) 



n=1 



where P}^ represents the associated Legendre polynomial. Figure [5] shows the modal de- 
flections /A for 2 < ri < 5 as a function of k. The simulation results are consistent with 
the measured deflect ions for all the shown m odes. In predicting the bubble migration, 
the wide-gap theory ( Magnaudet et al\ ( 20031) ) assumes that the mirror image primarily 
induces the deformation of the mode n = 2. The leading-order of the n = 2 deflection is 
/2 = K^/4 in the limit of k — ?■ 0. Considering the higher-order effect with respect to k, 



iMagnaudet et all (|2003t) derived fi^''^ = {n^ {1 + |k (1 + |k + g k^)}. For k <~ 0.7, 
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Figure 5. Modal deflection /4 ^\ which is expanded in associated Legendre polynomials P^, 
versus k(= (1 + £)~^). The dotted curve indicates the leading term of the analytical solution 
IjMagnaudet et a/1 (|2003f )) fn ~ \k^ , and the sohd curve the analytical solution with the 

higher-order contribu tions /A = jk ^ \1 -\ - |«; (! + §«:+ H i^^)}- The symbols indicate the 
experimental results (iTakemura et al\ (|2002l ')l. and the dashed curves the results obtained by 
the numerical simulation. 



such a higher-order k correction is responsible for the enhancement of the n — 2 de- 
flection from the leading-order one as seen in the better agreement with the measured 
and simulated deflections. However, the correction is not sufficient for the narrower gap 
K >^ 0.7, and thus the theoretical underestimation of the n = 2 deflection becomes 
more serious with n. Further, the theory does not cover the considerable increase in the 
higher-order n > 3 deflections with k as demonstrated by both the measurement and the 
simulation. 



The modal deflection /, 
into 



(Ca) 



is linked to the migration velocity and force as decomposed 



V, 



(Ca) 



B3 



^ y/Ca) 



B3,n' 



F, 



(Ca 
M 






(2.14) 



which are 



V, 



(Ca) 
S3,n 



Vbi 



p(Ca) 

Fdc 



r(Ca) 
RI,n 



d'x C 



(/^^''^^n (cos 6*) cos 



(2.15) 



(2.16) 



Figure |6l[a) shows the contribution of the modal deflection to the migration velocity for 
the modes n = 2 and rt = 3. For sraall k, the simulation result is consistent with the an- 
alytical solution (JMagnaudet et ali (|2Q03[ )'). which considers only the n = 2 deformation 
to cause the bubble migration. The inset shows that for small k, the contribution of the 
mode n = 2 is proportional to k^ , while that of the mode n = 3 to k^ , whose exponent 
is not trivially proven but may be predictable extending the regular perturbation to the 
higher-order. The difference in the exponent ensures that the relative contribution of the 
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Figure 6. Contribution of the modal deflection to the migration velocity Vg^J^, which is 
expanded in associated Legendre polynomials P^- (a) top panel: the linear-linear plot of the 
velocity of the modes n = 2 and n = 3 versus k{= (1 + s)~^). The solid line indicates the 

analytical solution (|Magnaudet et al\ pOOS)) V^'^'^^ /Vbi = 3k^{1 + 3k/2) /{40(1 + 3k/4)}. The 
symbols indicate the results obtained by the numerical simulation. The inset shows the same plot 
in a log-log scale, (b) bottom panel; The simulation results of the modal contribution normalized 
by the total velocity. 



mode n = 3 to n = 2 becomes more significant with k. It should be noticed that although 
we confirmed that the migration force contribution F|^ ^' of the mode n = 2 increases as 

K ^- 1 (i.e., e -^ 0), the velocity contribution Vgg 2 reduces as shown in figure[S{a). It is 

because the slope of the n = 2 migration force, — dlogFJ^ 2 /dloge, in a logarithmic plot 
is more gentle than that of the drag force, — dlogi^£)c/dloge — > 1 in the denominator of 
(|2.15p as £ ^- 0. Figure ini^b) shows the modal contribution V^g „ compared with the mi- 
gration velocity Vg^ . The contribution of the mode rt = 2 monotonically decreases with 
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K. The contribution of the mode n = 3 increases with k in the range of k <'^ 0.8, while 
decreases in the greater k. It is because the higher modal contributions for n > A become 
no longer disregarded. Moreover, the fact that all the shown modal contributions decay 
as K — > 1 indicates that the further higher-order contributions become considerable, and 
the regular perturbation approach with respect to k is no longer effective. 

2.5. Comparison with small deformation theory in the lubrication limit 

Figures [4l [5] and [6] imply that for small e, the bubble deformation is preferentially en- 
hanced within the narrow bubble-wall gap, and then its squeezing effect promotes the 
bubble migration. To shed more light on the role of the hydrodynamics in the gap, 
comparisons with a smal l deformation theory in the lubrication limit will be made. 
It should be noticed that iHodges et al\ ( 20041 ) have performed a lubrication study for 



a nearly spherical drop near a tilted plane in so-called 'slipping' regime, and derived 
the deformation-induced normal force. One can also access a relevant physical picture 
in theoretical studie s on t h e lift force on an elastic body induc e d by its deformatio n 
(|Sckimot o fc Leibled (|l993f) : ISkotheim fc Mahadevm] (|2004l[2005h : IUrzav et all (|2007l )V 
Following the spirit of the lubrication theory, we evaluate the migration velocity in the 
limit of £ — > 0. 

The basic equations for the lubrication analysis and the solutions are shown in Ap- 
pendix [Bj To be recalled and to be used for comparison with the simulation results, 
the preconditions and the perturbed quantities are detailed here. We prescribe the 
wall-parallel velocity Vbi, and employ standard lubrication assumption, i.e., e <^ 1. 
We also suppose small capillary number Ca ^ 1. As implied in (|2.2ip . the deflection 
is 0(Cae~^/^a), which has to be sufficiently smaller than the gap ea if the tilt an- 
gle of the near-wall interface from the plane wall is supposed to be small. Here we 
adopt an additional constraint S = e^'^/^Ca ^ 1. For inner coo rdinates (R, Z, ^) (her e 
e^/'^R = r and £~^Z = z a s illustrated in figure [ll[b), (see e.g. iGoldman et al\ (jl967[ ): 
O'Neifi fc StewartsonI (|l967l) )). a parabofic profile Z = H{R) = 1 + R'^/2 represents the 



interface within the inner region r r^ e^' ^, if the deformation is absent. Using e, wc write 
the velocity vector, the pressure, and the deflection f^'~^^> of the interface in an expansion 
form with respect to 6 



Ur{r, cj}, z) =U^°^{R, Z) cos (f) + 6Uy\R, Z) 

+ eUi:^\R, Z) cos (j) + SUy-'^\R, Z) cos2</- + OiS"^) + 0{e5) + 

u^{r,(t),z) =UJ^'{R,Z)sm(j) 

+ £U^^\R, Z) sin0 + SU^^-'^\r, Z) sin 20 + 0{d^) + Oie6) + . 

u^ (r, (/), z) =£1/2 fc/i") (i?, Z) cos + 5U^J^ {R, Z) 

+sUP{R, Z) cos(/> -f SUi^-^^R, Z) cos2(/) + 0(6^) + 0{e6) + .. 
p{r, 0, z) =£-3/2 fpio) (i^) cos + SP^f'> (R) 

+eP''^\R, Z) cose/} + SP^f'^HR, Z) cos2(j) + 0(6'^) + OieS) + ... 
/(^^) (r, <j>) = £-1/2 (f{R) cos (I) + 0{S) + 0{e) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 

whose scaling relations are suitable to all equations in Appendix |Bj It should be noted 
that the terms with the superscript (0), (1) or (/, 2) in (|2.17p - (|2.20p are proportional 
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to cos(f>, siiKJ), cos 2(/) or shi2(f), and thus provide no wall-normal force. Here we make 
clear the physical domain of validity of the condition S <^ 1 when the rising velocity 
Vbi, which is needed t o eva luate Ca(= hVbi/'j), is unknown. Following the analysis by 
O'Neill fc StewartsonI ( 19671 ). we write the drag force Fd acting on the spherical bubble 



translating parallel to the wall in a form Fd — 6tt ^aVBi{A\og e + B) , where A and B are 
independent of e. Considering the free-slip boundary condition on the bubble surface, we 
can analytically find A = —1/5. From our numerical data, we approximately estimate 
B = 0.6. Hence, from the force balance Fd = Anpa'^g/S, for e ^ 1, we evaluate the rising 
velocity 

1 lOpa^g 

""' (-log£ + 3) 9/i ' 
and obtain the following relation for the bubble radius a to satisfy the condition (5 <C 1. 



To assure us of the appearance of the lubrication effect, figure[7]shows the profiles of the 
pressure p(*')(= 3R/{5H'^)) and defiection -F(= 31ogiJ/(5i?)) in the lubrication limit 
(see (jB 11[) and ()B 12p . respectively), which are compared with the simulation results of 
the scaled interfacial pressure e'^/'^P/ cos(f> and the scaled deflection £i/2j(Ca) j^gg,!- the 
wall as a function of £~^/^0. For sufficiently small e, the simulation data of the scaled 
pressure collapse onto the curve (|B 11[) as shown in figure [Tl^a). Therefore, the scaled 
defiection profile reasonably approaches the lubrication solution (|B 12p with decreasing 
the bubble- wall gap as shown in figure [Tljb). As plotted as the dashed-dotted curve in 
figure m the defiection based on (|2.2ip and (|B 12|) is consistent with the measured and 
simulated deflections in the wall neighbour [6 ^ 0). Hence, the dis crepancy between 



the de flections of the narrow-gap experiment and the wide-gap theory (jMagnaudet et 
( 20031 )) is attributa ble to the lubrication effect. 



From lBartI ( 19681 ). the drag force acting on the bubble translating perpendicular to the 



plane wall is Fdc — > 37r^/(2e) as e ^^ 0. Substituting this relation and (|B20|) (i.e., Fm 
in the lubrication limit) into \2.9\ , we obtain the asymptotic solution of the migration 
velocity 

To make comparisons with the asymptotic solutions, figures [2] and [8] show the scaled 
migration velocity as functions of the inverse distance k(= (1 -|-e)~^) and the clearance 
parameter e, respectively. The inset of figure [8] shows the scaled force (jB 20[) . The simula- 
tion results are consistent with two asymptotic behaviors based on the lubrication theory 
for e <^ 1 as well as the mirror image technique for e ^ 1 (i.e., n <^ 1). As seen from 
figure m the theories provide the different exponents of the migration velocity scaling 
with respect to e, namely, Vsz/yBi on Ca e^^ for e :s> I and Vb^/Vbi oc Ca e^^ for 
e <C 1. 

The puzzling finding in iTakemura et all ( 20091 ) that the theory of lMagnaudet et 



(|2003( ) accurately predicts the deformation but fails to predict quantitatively the deformation- 
induced migration velocity is explained from the fact that the ratio of the simulation 
result of Vb3 to the theoretical prediction is more sensitively dependent upon e than 
that of /j when the lubrication effect becomes relevant. For instance, for e = 0.4, 

"(Ca) 

£ = 0.2 and e = 0.1, the simulation-to-theory ratios of /j are respectively 1.004, 1.1 
and 1.3 (figure [5|), while those of Vsa are respectively 2.5, 4.3 and 7.3 (figure [2]). The 
lubrication effect is likely to compensate the large discrepancy between the migration 
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Figure 7. Profiles of the pressure ((a) top panel) and the deflection ((b) bottom panel) in 
the inner expansion scalings of the lubrication theory. The solid cu rves indicate the predictions 
HB ll|l and (|B 12^ by the lubrication approach (jHodges et al\ (|2004f l'). The symbols indicate the 
results obtained by the numerical simulation. For the plot, the simulation data of the interfacial 
pressure P and deflection f^'^'^' for the bubble translating parallel to the wall are scaled as 
pm ^ e-3/2p/ cos and F = e'^''^ f^'^''\ and the angular coordinate as 7? = e'^^'^O. 



velocities of the experiment and the wide-gap theory revealed in lTakemura et aU (|2009| ). 
However, although the quantitative agreement between the interfacial deflections of the 
experiment and the simulation is shown in figures |4] and [S] the migration velocity in the 
experiment is still considerably higher than the simulated one. Its cause is not clear at 
this moment, and further joint researches among theory, numerics and experiment are 
needed to resolve this discrepancy problem. 
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Figure 8. The lateral migration velocity versus the clearance parameter e. The circles indicate 
the results obtained by the numerical simulation. The solid, and dashe d-dotted curves correspon d 
to the analytical solution Vbs = SCaVsiK^ (1 + 3k/2)/{40(1 + 3k/4)} (iMaenaudet et al\ (|2003f) '). 



and Vb3 = ^CaVsie ^ derived by means of the lubrication ap proach (iHodges et al. ( 2004 )). 



respectively. The dashed curves indicate the experimental results (JTakemura et all (|2002l )V The 
inset shows the lateral force (|B 20|l . 



3. Conclusion 

We numerically and theoretically investigated deformation-induced lateral migration 
of a bubble slowly rising near a vertical plane wall in a stagnant liquid. We focused 
on a situation with a short clearance c between the bubble inter face and the wall. We 
demonstrated that the wide-gap theory ( Magnaudet et all ( 20031 )). which considers the 
n = 2 deformation mode, describes the deformation-induced lift force as long as the 
bubble-wall gap is sufficiently wide {a/{a -|- c) ^ 1, here a is the bubble radius). For 
the narrow-gap case with the clearance parameter £(= c/a) smaller than unity, we found 
that the higher-order n > 3 d eformation mode s crucially enhance the migration velocity, 
and the lubrication effect (jHodges et all ( 20041 )) appears to induce the migration velocity, 
which scales asymptotically like Vbs — > Ca e^^VBi/25 as e — > 0. This contrasts with the 
case of the inertia-driven migration, to which the wide-gap theory demonstrated a robust 
applicability in prediction over a wide range of e. 

The present simulation consistently served as bridge between the wide- and narrow-gap 
theories (see figure [8|) as long as the bubble deformation is assumed to be infinitesimal. 
However, in spite of the qualitative success of the simulation in revealing the narrow-gap 
repulsive force, the deformation-induced migration velocity in the experiment is con- 
siderably higher by a factor of about 3 than the simulated one, as shown in figure [H 
The experiment may inevitably involve unknown factors such as unsteadiness, imper- 
fection from the infinite plate-fluid system, and measurement uncertainty, which cannot 
be captured by the simulation. Nevertheless, we have not expected such a large dis- 
crepancy because firstly the quantitative agreement between the interfacial deficctions 
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►^ 

Figure 9. Computational grid, (a) left panel: The bipolar coordinates. ^ = and ^ = a 
represent the plane wall and the bubble surface, respectively, (b) right panel: Definition points 
of the velocity components and the pressures on the staggered grid in a computational space. 



of the experiment and the simulation has been confirmed in figures |4] and [S] secondly 
the c onsistent inertia-driven migration velocity of a r igid sphere with the av ailable the- 
ories (IVasseur fc Coxl (J1976I ): iMagnaudet et al\ pOOSh ) has been obtained bv iTakemural 
( 2004^ using the same experimental setup, and thirdly considerable uncertainty seems 
not to be introduced into such a simple system as illustrated in figure [TJ Further joint 
researches among theory, numerics and experiment are needed to resolve this problem. 
As a possible factor causing the inconsistency, we raise a difference between the bubble 
deformation levels of the experiment and the present analysis. As stated in §2.51 the 
infinitesimal deformation assumption in the lubrication limit is justified only for the case 
that 5{~ e~^/^Ca) ^ 1. Beyond this limitation, unexplored hydrodynamic ingredients 
possibly become important on the bubble migration. For the experimental data shown in 
figure[2]and figure[8l the maximum value of 5 is 0.74, which is less than but comparable to 
unity. Therefore, the bubble deformation is finite rather than infinitesimal, and is likely 
to induce the higher-order force, which is possibly comparable to or stronger than the 
leading migration force evaluated with the infinitesimal deformation theories. From the 
theoretical viewpoint, a tiny bubble experiment, which results in a tiny capillary number 
and thus a tiny deformation, is favorable for comparative study. However, such an exper- 
iment has often resulted in an undetectably low migration velocity and made an accurate 
measurement difficult. T o overcome such a dilemm a , the h ighly accurate boundary ele- 
ment computations (e.g. IWang &: DimitrakopoulosI (|2006l ) : IPimitrakopoulosI (|2007l )) for 
various deformation levels would be helpful to complement the infinitesimal deformation 
theories. 



We thank Shu Takagi for fruitful discussion. 



Appendix A. Finite-difference descriptions of the basic equation set 
in bipolar coordinates 

We describe the basic equation set in the bipolar coordinates (see e.g. lO'Neilll ( 19641 ): 
Happel fc Brenneii (|l973[ ) Appendix A-19) as illustrated in figure [9lja). The coordinates 
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(r, z) in figure [1] (b) are 
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fcsinyy 



k sinh ^ 



(Al) 



where 2? = cosh £_ — cos 77 and k = a/£(F+2) . The bubble surface is located at ^ = a = 
log(l + e + k) as shown in figure [Hla). The gradient of a scalar function q is written as 



Vg = 



Gf 



£0 5g 



(A 2) 



where e represents the unit vector, and its subscript the corresponding component, h 
denotes the scale factor defined by e.g. h^ = {{dxi/d^)^ +{8x2/8^)^ +{8x3/8^)^}^^^ 
(see e.g. lBatcheloiJ ( 19671) . Appendix 2). Each component is explicitly given by 



K 



k 



(A3) 



The vector and pressure field ([/, P) for the bubbl e translating parallel to the w all is 
written using the quantities with hat in a form fe.g. lSugivama fc Sbragaglial ( 20081 )) 



U = 



{{eA 



S^-jU j-j 



cost 



e^f/^ sin (/) — ei ^ , P = P cos ( 



(A 4) 



of which the Fourier expansion reduces the three-dimensional problem to the two-dimensional 
one. The vector and pressure field (m,p) for the bubble translating perpendicular to the 
wall is 

M = e^uj -I- e^u^ + 63, _p = p. (A 5) 

It should be noticed that as opposed to the reference frames in ij2.1[ we take those for 
(t/^, Un, U^) and (u^, u,,) viewed from the plane wall for convenience of the simulations. 
Consequently, U^ = Uri = U^ = u^ = u,-i = on the plane wall and on the boundary 
sufficiently far from the bubble. 

We follow a conventional staggered grid arrangement (JHarlow fc WelchI (|l965l) ). where 
the velocity component is located on the corresponding cell interface, and the pressure 
at the cell centre, as shown in figure |9l[b). The basic equations are discretized by the 
second-order finite-difference scheme. To numerically guarantee the mass and momentum 
conservations, and to accurately conduct the numerical integration in computing the drag 
and migration forces, we use the exact values of the grid width, the cell intcrfacial area, 
and the control volume in the bipolar coordinates. To this end, for the integral of the 
scale factors 

9d^,v)^ dCh^{C,v), 9i4>i^,v)^ <iC hi^^v)h4,{^,v), 



9ivi^^v) = / d^ / dr? h^i(,fi)hr,{^,fj), g{C,v) = / dC / df/ h{^,fj), 



a J 7T 



a J TV 



we use the exact expressions 



2k ( _i (V + C\ ^ _i (V^+Ca. 
9^ = — — I tan I — - — I — tan ' 



sin r]\ \ S 



2k _, fV + C 

9ri ^ r-r-: tan 



Oq 



sinh ^ 



S 



(A 6) 

(A 7) 



g^n =fc- 
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2\ f I 1 \ ,____-! fV + C' 



\sinh ^ sin rj 
(I 1 



tan 



S 



tan 



Vsinh a sin 77 

gr,0 = fc^ ( -7^ 



i(Va+Ca\ f C + 2 Ca + 2 
2S 2Sa 



1 1 



9i4> = k 1-— 2— (tan \—g— 



sin ?7 



V cosh^ + 1 

1 / 2?Q, + Cq \ \ sinh ^ sinh a 
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(A 8) 

(A 9) 



tan 



5a 



T) sin 7] I?Q sin 77 



sin 77 



tan I — :; — I — tan | | | cos 77 -\ 



where 



S J \ Sa J J 22? 22?, 

tanh'' - - tanh'^ - H tanh - - tanh - 

12 V 2 2y4V2 2 

h = h^hrihif,, C = cosh^cos77 — 1, 5 = sinh ^ sin 77, 

"Dq = cosh a — cos 77, Cq = cosh a cost; — 1, Sa ~ sinh a sin 77. 



(A 10) 
(All) 

(A 12) 



We introduce the finite-difference operators Si and Sj , of which the indices i and j cor- 
respond to discretized coordinates along the respective directions ^ and rj, such as 



<5j(5j((7)|,j 



9«+i,j+i -'?»+i,j-i -9^-ij+i +'?»-ij-i- 



(A 13) 



Using these operators, we write the divergence of the velocity vector U in (|2.ip as 



V-U 



i,j \ coscp / SiSj{g)\ij 

(A 14) 
and the components of the Stokes equation (|2.2p for the bubble translating parallel to 
the plane wall as 



= 
= 
= 



(5j(5?)|i+ij Sj{gr,^)\i+ij 



(P-V^)l,,, , -(5,((5j(g,,)i7^)|,,, +(5J((5,(54)i75)|,,, 



l'»J 



Si6j{g^,^)\i^j 



(A 15) 



where J7 denotes the vorticity, of which each component is 

S3{rU4,)\ij+i + UvSj{g^)\,j+i 



n 



iKj+i 



n 



^\i+hJ 



^0lj+ij + i 






(A 16) 
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Replacing U^, Un and P respectively by u^, u,, and p with C/<^ = fij = $7,, = 0, we 
readily obtain the governing equation for the field {u^^UriTP)- We write the kinematic 
and free-slip conditions (|2.4p as 

^?|jv„ + ij == er^|jv„+ij, ucU^+i J = -e^ • e^l^^+i J, (A 17) 

where the index Na + ^ denotes the node at the bubble surface ^ = a, 



5 



-,(«) 



(•9eV) . „ _ '5,(r2) 



^^■^^ = ^;m' ''-'^^-w^y ^^''^ 



m r, / k^Ssmfi\ ,o 1 ,(V + C\ Csin?7 , , _, 



■^4') 



(5»(^j(ffC'j) 



'5.(9c) '5j(C/c)-^^'5j(5c)f^C + 5,{9r,) 5,{U,,) ~ 5,5j{g^)U^ , (A21) 



S?0 = j-^^ |r^5,(C/^) - 5,(r)C/^ - 5^{g^)U^\ , (A 22) 

o-e-? = ^.^ ./ N {*»(ffc)^'5j(tic) - '5,5^(3^)^5^ +5j{grSHur,) ~ 5^5J{gr,)un'^ , (A23) 
and the overline stands for the interpolation such as 

We write the area integral on the bubble surface in a summation form 

(tdx^ q= d(f) d?y h,jh^ q\^=a = X! ( ^3^3,i<p) / d</) g ) , (A 25) 

"'Sb Jo Jq j^-^ \ Jo / w<,+i J 

where Nj is the number of grid points in the 77-direction. The drag force Fdq in (|2.10p 
is given by 

Fdc = 27r ^ Uj ig^4,)e^ ■ e^a^A , (A 26) 



'-^ + i.' 



where 



For the deflection / ~ /''-''^^/(acosc/)), the Laplace law (|2.3I) is expressed as 

ai(]/|,+i + al{,Vl,_i - a({)/|, = ^^-^'l,, (A28) 
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where 






Af) 



rjv +i 1 + i 



Sji9ri) U^+i j+il5j (5^0)1 Af^ + j 



-2 J+2 

JVq-T 2 ij 2 






P,J 



{er■e^)^^_^_lJ 



(A 29) 



,iV. 






(A 30) 



^ii = -P 



2Sj{g^^) S,{U^) , 2rJ<5,<5,(g5)[/^ 



6,6, (g) 



The deformation-induced lateral force FJ^j ^' in (|2.7p is given by 






6,6 jig) 



(A 31) 



F; 



(Ca) _ 



M 



-E 



-»^i( 



^j{9v4,)^ii h{Ui - e r-e^)f r6,6j{g ^)a^^ (C/^ - e^-e^) / 



-^iCg?) 



-^iCffs) 



+ Sj{g^4,)^a' — ^- / - ra^i' (Ur, - e^^) ^j(/) 



Hi - Hn 



Si (90 



S^i90'' 



+ r{urj + e^-e,,) S,,,, 6j{f) - 



yUrj \ 6^ ■ ^7} ) 



5.j{gn^)T.,-iA f 



+ r(u,, + e,-e^)5(/) 5j{f) 



iV„ + i 



/ 



(A 32) 



where 



An) 

9^r, 



-r ^7] 



de 



SjigSl 

Si{9i4>) ' 
1 






j,^ fc^CsinT? _ , 2 j 

2?^ I sin^ 77 



tan 



1 /2? + C 



5 

Ca sin 7] 

25 25a ) V22?2sinhC 2P2 sinha 



tan 



1 ( ^a I ^Q 



C + 2 Ca + 2\ / C sin 77 



(A 33) 



(A 34) 



^nn 



-P 

1 



26,ig^^y 6,iU,j) , 2r6,6,ig,)U^ 



5,6 j{g) 



6^6 jig) 



-"<i4> 



Sj{9v<l>) 



r'6,{U^)-5,{r)U^ - 6,{g^)U, 



(A 35) 
(A 36) 
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Appendix B. Small deformation theory in the lubrication limit 

The governing equations for U^"\ [/'"', ui"\ p("), f/^', uP and P^^^ in ^JT^-^M^ 



are written as 



1 d{Rm 

R as" 



(0)x 



u. 



(0) 



dU: 



R 



dZ ~ R OR ^ dZ 



(Bl) 



,2r>(0) 



dR 

dpU) 



dZ^ 



p(0) 



dR. dZ^ dZ 

with the no-shp boundary condition on the plane wall 

M = at Z = 0, 



dZ^ 
= 0, 



apw 
az 



(B2) 



(B3) 



and the free-slip and kinematic boundary conditions on the bubble surface. The deformed 
interface is located on the curve where H — 5F cos (j) — Z = holds. From this relation, the 
normal unit vector n pointing outwards the liquid on the bubble surface is approximated 

by 



n{R, 



1 /2 / dF 

ErE ' \ R— 0-— -C0S( 

\ dR 



e^e^''^5— sin( 
R 



(B4) 



Applying the Taylor expansion to a function g = q^^^^ + 5 q^^^^ + . . . in terms of the deflection 
around the undeformed interface, one obtains a relation on the deformed interface Z = 
H -SF cos 6 



9 1 interface =9lz=_ff +^ { - F COS 



dqiO) 



dZ 



1\z=H 



Z=H 



(B5) 



Taking (|B4[) and (|B5|) into account, one writes the kinematic condition on the bubble 
surface as 






FdU. 



(0) 



2 dZ 




and the free-slip boundary condition as 



dU^'^ dU'f^ dU^r 



rif) 



(0) 



dZ 



dZ 



dZ 



Fd^U, 



2 az2 



at Z = H. 



(B6) 



(B7) 



(Bi 



The vertical drag force acting on the bubble, which is involved in the second term in 
the right-hand-s i de of (I2.3p. is of order loge as determ ined for a motion of a rigid sphere 
(JGoldman et al\ (|l967l ) ; lO'Neill fc StewartsonI (|l967l) ) by means of the matched asymp- 
totic expansion technique and the normal stress on the bubble surface is dominated 
by the pressure p = 0{e~^/^) as compared with drUr = 0(e~"^/^), Ur/r = 0{e~^^'^), 
u^/r = 0(e~^/^) and d^Uz = 0{e~^^'^), which are related to the viscous stresses. Hence, 
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the Laplace law (|2.3|) is simplified into 

dRyn dR j ' ^ ' 

with no singularit y conditions f = a t R — and F — > as i? — > oo. 

As obtained bv lHodges et all ( 2004 ). the leading-order and perturbed solutions are 



~ (0) _ (6 - 9fl^) p^ 2Z\ fA0)^_3_fZ^_2Z 
20H \H^ Hr '*' 10li72 H 



m 



{AR - R^)Z'^ {-42R + iR^)Z^ 



(BIO) 



5if4 2{)H^ 



P'»'^||. (BU, 

F = ^. B12 

5i? ^ '^ 

Uy\R,Z) ^A2{R)Z^ +Ai{R)Z, 

m/HR,z) = -^^-^^^z^-^^-^^^z^ ^^''^ 

^ ^ ' '^ 3R dR 2R dR 

rR _ 
P^f\R)^ / dR2A2{R), (B14) 

9(14-i.2)logi^ (B15) 

Substituting (jB 16p into (jB 14p estimates the asymptotic order of the perturbed pres- 
sure 

which is 0(Calog£ (e^l'^K)^^) in the overlapping region R ^ e"^/^, and thus to be 
matched with the pressure of 0(Caloge r~^) in the outer region (|0'Neill fc Stewartson 



where 



( 19671 )). This outer pressure may contribute to the lateral force of O(Caloge), which is 



la rger than tha t for th e wide-gap case, corresponding to 0(Ca). Nevertheless, as discussed 



in lUrzav et al\ (J2007r ). the rapid decay of the perturbed pressure for i? :§> 1 indicates 
that the contribution of the outer pressure to the lateral force is negligibly smaller than 
that of the inner pressure. Therefore, to evaluate the leading-order migration force, one 
does not have to solve the outer problem. (In fact, the order of the lateral force (jB 20p 
evaluated only in the inner region is confirmed to be 0{Ca£~^) and larger than the outer 
contribution 0(Calog£).) The leading pressure e^'^/^P^"^ cos with no deformation is 
locally dominant but does not contribute to the lateral force due to its azimuthal cosine 
dependence. The viscous stress contribution to the lateral force is 0{e) smaller than the 
pressure in the inner region. The deformation-induced lateral force Fm to cancel the 
migration velocity and to maintain the wall-parallel motion is expressed as 



/ d2a;e-3/2<5p(/). (B17) 



Fm 

^contact 
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The surface integral for a function q on the contact side is taken from the axis i? = to 
the overlapping region, i.e., 

f d'^xq = e d(t> dRRq{R,(t>), (B 18) 

^contact Jo JO 

where *H = 0(e^^/^). For y{^ 1, one obtains an asymptotic relation 

pfyK /-SH /-oo 

/ di? i? P(/) (i?) = - / di? i?2 ^2 (i?) - 9^2 / d^ ^2 (i?) 



3 , 61og(fHV2) , 1 , oi^-Ho,^l 



(B19) 



100 25«n4 59^4 

of wh ich only the first term does not vanish as ?l — ?► oo. As obtained by iHodges et al 



(200J), consequently, the asymptotic solution of the deformation-induced lateral force in 
the lubrication limit is 



r.2Ti- /.JR 



Fm -^ lim e <l(j) UrR e'^/^SP'^^'^R) 

3^^°° Jo Jo (B20) 

37rCa _, 
— - as e — )► 0. 
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